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Abstract 

It is known that in the zeta function regularization and in the Fujikawa 
method chiral anomaly is defined through a coefficient in the heat kernel 
expansion for the Dirac operator. In this paper we apply the heat kernel 
(D , methods to calculate boundary contributions to the chiral anomaly for 



' local (bag) boundary conditions. As a by-product some new results on the 

heat trace asymptotics are also obtained. 
^ : PACS: 11.15.-q, 11.30.-j, 02.40.-k 
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1 Introduction 

Chiral anomaly, which was discovered more than 30 years ago fT^, still plays an 
important role in particle physics. On smooth manifolds without boundaries 
many successful approaches to the anomalies exist^. Modern developments in 
theoretical physics require anomaly calculations on branes and domain walls (see, 
e.g., jE]). Taking into account also more traditional applications as the bag model 
of hadrons [7] , we see that understanding the chiral anomaly in the presence of 
boundaries or singularities is an important task. 

*eniail: inaraval@mail.ru, root@VM1485.spb.edu 

^Also at V. A. Fock Insitute of Physics, St. Petersburg University, Russia; email: 
Dmitri . Vassilevich@itp . uni-leipzig. de 

^The reader may consult, for example, ref. pj or the original papers 
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In the case of non-trivial background fields, and especially in the presence 
of boundaries or singularities, the heat kernel technique seems to be the most 
adequate one for analysing the one- loop effects (see JHj for a recent review) . The 
heat kernel approach to the anomalies is essentially equivalent to the Fujikawa 
approach 4j and to the calculations based on the finite-mode regularization [Oj, 
but it can be more easily extended to complicated geometries. 

In the present paper we consider an euclidean version of the bag boundary 
conditions [7j (see also JO])- Although these are the most simple and physically 
most natural boundary conditions, chiral anomaly for these conditions has not 
been calculated so far^. This can be explained by the fact that the calculation is 
indeed rather involved. We use the zeta function regularization and our approach 
to the anomaly follows closely the paper 11 . Mathematical foundations for 
analysing spectral geometry of the Dirac operator with local boundary conditions 
were developed in [12j. A simple overview of spectral properties of Dirac operator 
on manifolds with boundary can be found in [TB] . 

To calculate the anomaly we employ the following strategy. First we relate the 
anomaly to a heat kernel coefficient for the square of the Dirac operator. Then 
we generalise the problem and consider the heat kernel for an arbitrary operator 
of the Laplace type with mixed boundary conditions and with a matrix valued 
smearing function. This generalisation allows us to use the heavy machinery of the 
heat kernel expansion for such operators. In particular, important information 
follows from two simple examples of boundary value problems in one and two 
dimensions. In this way we are able to calculate first five heat kernel coefficients. 
Chiral anomaly is then obtained simply by substituting explicit expressions for 
the connection, the potential, and other quantities in terms of the background 
vector and axial vector fields. 

In the next section we remind some basic formulae regarding the zeta function 
regularization, chiral anomaly, and local (bag) boundary conditions. Sec. 3 is 
devoted to the relationship between Dirac and Laplace operators. The heat 
kernel coefficients are calculated in sec. 4, which is the main technical part of this 
paper. In sec. 5 we return to the Dirac operator and finally calculate boundary 
contributions to chiral anomaly in two and four dimensions. Our results and 
their possible extensions are discussed in sec. 6. Appendix ^ contains the heat 
kernel coefficients for mixed boundary conditions with a scalar smearing function. 
Appendix ^ gives details of the particular case calculations used in sec. 4. 

^There exist calculations of global chiral anomaly (i.e. for a position-independent chiral 
transformation parameter) for some other types of boundary operators, and also calculation of 
the scale anomaly for bag boundary conditions. Literature on these topics is rather large, so 
that for uniformity we quote nobody here. A literature survey can be found in [S]. 
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2 Boundary conditions and chiral anomaly 



Let us consider an n-dimensional Riemannian manifold M with boundary DM. 
Dirac 7-matrices satisfy the Chfford commutation relation 

{7^ 7^ = -2^?'^^- (1) 

The 7-matrices defined in this way are anti-hermitian, 7^^ = —7'^. We also 
need the chirality matrix which will be denoted 7^ independently of the dimension 
of M. As usual, 7^7^^ = —7^7^ and 7^''' = 7^. We assume that n is even. We fix 
the sign of 7^ by choosing 

•n(n+l)/2 

7' = ^e'^'-'l.lu.-.lp. (2) 

Consider the Dirac operator 

D = r {d, + V, + zA^75 - \bp, iMr^^ (3) 

in external vector and axial vector fields. We suppose that and are 
anti-hermitian matrices in the space of some representation of the gauge group, 
crjf^^ is the spin-connection^. 

The Dirac operator transforms covariantly under infinitesimal local gauge 
t r ansf ormat ions : 

5xA^ = [A^, A] 
SxV^ = d,X + [V,,X] 

D^D + [D,\] (4) 

and under infinitesimal local chiral transformations: 

6^Af, = df,ip+ [Vf,,ip], 
5^V^ = -[A^,(p], 

D ^D + i{D,^''^]. (5) 

The parameters A and are anti-hermitian matrices. 

We adopt the zeta-function regularization*^ and write the effective action for 
the Dirac fermions as 

ly = - In det 5 = In det = ^('(0) + \ H^^)m , (6) 

•^The spin-connection must be included even on a flat manifold if the coordinates are not 
Cartesian. 

''Our approach to the effective action for fermions is close to that in the papers jlll I14j . 
We refer to these works for a more detailed derivation of chiral anomaly (with somewhat more 
accurate treatment of the zero mode problem). 
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where 

C(.) = Tr(5-2.) ^ (7) 

prime denotes differentiation with respect to s, and Tr is the functional trace. 

It is easy to show that the effective action (jHj) is gauge invariant, 6\W = 0, 
and that the variation of W under an infinitesimal chiral transformation reads 

A := ~6^W = -2zTr(7V5-2^)|,=o • (8) 

Let us define an integrated heat kernel for a second order elliptic partial 
differential operator L by the equation: 

KiQ,L,t) :=Tr(Qexp(-tL)) , (9) 

where Q{x) is a matrix valued function. For the boundary conditions we consider 
in this paper (see eq. below) there exists an asymptotic expansion [19] as 
t 0: 

oo 

K{Q, L, t) ~ ^ ak{Q, _ ^^^^ 

k=0 

The heat kernel is related to the zeta function by the Mellin transformation: 

/■oo 

Tr(7V5"'') = r(s)~M dtf~'K{^^ip,D'',t). (11) 

Jo 

In particular^, 

^ = -2m„(7V,52). (12) 
The same expression for the anomaly follows also from the Fujikawa approach 

We shall need some basic notions from differential geometry. Let R^upa be 
the Riemann tensor, and let i?^,^ = fj.ua be the Ricci tensor. With our sign 
convention the scalar curvature R = R'^^ is +2 on the unit sphere S*^. Curvature 
does not play any important role in our calculations. However, we shall see below 
that curved space offers no complications in our approach compared to the fiat 
case. 

If the manifold M has a boundary, boundary conditions should be imposed 
on the spinor field ip. We need several basic definitions regarding differential 
geometry of manifolds with boundary. Let {e-, }, j = 1, . . . , n be a local orthonor- 
mal frame for the tangent space to the manifold and let on the boundary e„ be 
an inward pointing normal vector. Then {Ca}, a = 1, . . . , n — 1 can be identified 
with a local orthonormal frame for the tangent space to the boundary. The frame 
{cj} will be used to transform curved (world) indices /x, z/, . . . , a to "fiat" indices 
and back. For example, for a vector this transformation reads: vj = e^f^, 

rather formal way to derive this equation consists in integration of the asymptotic ex- 
pansion l|10() . 
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Va = e^f^, Vn = ef^Vfj,. In Euclidean space there is no distinction between flat 
upper and lower indices. 

The extrinsic curvature is defined by the equation 

Lab = ^ab 5 (13) 

where F is the Christoffel symbol. For example, on the unit sphere S"""^ which 
bounds the unit ball in M" the extrinsic curvature is Lab = ^ab- 
We impose local^ boundary conditions: 

n_v^|aM = o, n_ = ^(i-757„) , (14) 

which are nothing else than a Euclidean version of the MIT bag boundary con- 
ditions [7]. For these boundary conditions III = n_, and the normal component 
of the fermion current V'Wn^ vanishes on the boundary. 

An important comment on chiral transformations of the boundary conditions 
()14|) is in order. Finite version of the infinitesimal transformation (0) reads: 

D^ = e'^^^We'"^^' . (15) 

This relation yields the following transformation law for the boundary projector: 

n_ ^ n^:^^ = e-'^^m_e''^^^ = - (l - TsTne^^*^^') , (16) 
so that the boundary condition 

nL^l^laM = (17) 

remains consistent with (fT^ and (fT3j) . Eq. (fT7|) represents an Euclidean version 
[T^ of chiral bag boundary conditions ^7]. The boundary conditions (fTTj) are 
considerably more complicated than (jl4p . Even such fundamental property of 
(|17|) as the strong ellipticity (which ensures, for example, existence of only simple 
poles of the zeta function) has been established only recently ^HI- Fortunately, as 
we stay at the level of linear perturbations, the condition (fT^ is enough for our 
purposes. However, already the Wess-Zumino consistency conditions which 
imply two consequent chiral transformations require a more general setting of the 
chiral bag (fTTj). 



3 Dirac and Laplace operators 

To calculate the anomaly (fT^ it is convenient to consider a more general problem 
of calculation of the coefficients ak{Q, L) (see eq. fjlO|) ) for general matrix valued 

^Locality means that the projector n_ acts at each point of the boundary independently. 
An example of non-local boundary conditions can be found in |15| . 
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function Q and general operator L of Laplace type. Any operator of Laplace type 
can be expanded locally as 

L = -{g^'^d,d, + a''d, + h), (18) 

where a and h are some matrix valued functions. One can always introduce a 
connection uo^ and another matrix valued function E so that L takes the form: 

L = -{g'''V,V, + E) (19) 

Here is a sum of covariant Riemannian derivative with respect to metric g^^^ 
and connection u)^. One can, of course, express E and a; in terms of a^, h and 

^M = ^^7M.K + ^?'^r;j, (20) 
E = h-g^'''{d,u^ + u^uj,-UpT''^,) (21) 

For the future use we introduce also the field strength for uj: 

^iiu = d^ujy - duUJf, + [ujf,, Uu] . (22) 

The connection cj^ will be also used to construct covariant derivatives. It will 
be convenient to work with fiat indices (denoted by Latin letters) which we have 
introduced in the previous section. The subscript ]i . . .jk will be used to denote 
repeated covariant derivatives with the connection u and the Christoffel connec- 
tion on M. The subscript : a . . .be will denote repeated covariant derivatives 
containing u and the Christoffel connection on the boundary. Difference between 
these two covariant derivatives is measured by the extrinsic curvature ()1H|1 . For 
example, E-ab = E,ah - LabE-n- 

Let us now turn to boundary conditions. We assume given two complemen- 
tary projectors n±, n_ + n_|_ = / and define mixed boundary conditions by the 
relations 

n_7/'|aM = 0, (V„ + 5)n+^|aA/ = 0, (23) 

where S" is a matrix valued function on the boundary. In other words, the compo- 
nents n_'?/' satisfy Dirichlet boundary conditions, and Tlj^ip satisfy Robin (modi- 
fied Neumann) ones. 

It is convenient to define 

x = n+-n„. (24) 

Now we have to calculate the geometric quantities introduced above for generic 
L in the particular case L = D^. After lengthy but straightforward calculation 
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one obtains from Q, (PU)) - (see also [22] for the abelian case): 



u^, = V,- ^[7„ lAA^l, - iMr^, (25) 

E = -\ri''V^. + {n- 3)7^"A^A, - A^A^ --^R + iD^A'^^,, (26) 
= t-f%D^A^)j^j5 - t^''{D^A^)-f^^5 + iA^^-f5 - [A^, A^] + V^^ 




with the notations V^y = (9^K - d^Vf, + [V^, K], Af,^ = D^^A^ - D^Af,, D^A^ = 
dfj_Ay — T^^j^Ap + Note, that the gauge covariant derivative differs 

from defined above. 

Since D is a first order differential operator it was enough to fix the boundary 
conditions ()14|) on a half of the components. To proceed with a second order 
operator L = we need boundary conditions on the remaining components as 
well. They are defined by the consistency condition: 



which is equivalent to the second (Robin) boundary condition in (j23|) with 



4 Asymptotic expansion of the heat kernel 
4.1 General strategy 

In this section we study the short t asymptotics (jlUj) for an arbitrary operator L 
of Laplace type. Two particular cases of the expansion (jl(J|) are known. The heat 
kernel coefficients ak, k = 0,1, 2, 3, 4, 5 for a scalar Q = fl (where / is a function 
and / is the unit operator) are presented in Appendix ^ The case of arbitrary 
Q but pure Dirichlet or Neumann boundary conditions (i.e. when either 11+ or 
n_ is zero) was studied in |2S]- Here we need a combination of these two cases. 
Namely, we are interested in mixed boundary conditions and a matrix valued Q. 

According to the general theory \XQ[ the coefficients a^iQ, L) are locally com- 
putable. This means that each ak{Q,L) can be represented as a sum of volume 
and boundary integrals of local invariants constructed from Q, Q, E, the curva- 
ture tensor, and their derivatives. Boundary invariants may also include S, Lab 
and X- Total mass dimension of such invariants should be k for the volume terms 
and k — 1 for the boundary ones. 



n_D^|aM = 0, 



(28) 



S 




aa 



(29) 
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The following property of the heat kernel coefficients ^H] will be useful in the 
calculations. Let us define a shifted operator 



L, = L-eQ. (30) 

Then 

Tr(exp(-tL,)) = tTr(g exp(-tL)) . (31) 



d 
de 



e=0 

By expanding both sides of this equation in a power series of t one obtains: 



d_ 

de 



ak+2{l,L^) = ak{Q,L) . (32) 



e=0 



All geometric quantities (metric, effective connection, boundary conditions, etc.) 
corresponding to are the same as for the unperturbed operator L except for 
the potential E which receives a shift, 

E, = E + eQ. (33) 

Therefore, variation of the heat kernel coefficients w.r.t. e is equivalent to the 
variation of E. Equation together with the heat kernel coefficients for scalar 
smearing function presented in Appendix |X] allow to calculate the coefficients 
ak{Q,L) for k = 0,1,2,3: 



ao(Q,i:) = (47r)-'^/M rfna;y^tr(g). (34) 
Jm 

a^{Q, L) = ^(47r)-("~i)/^ / d^-^x^ti (xQ). (35) 
4 JdM 

a2iQ, L) = ^(47r)-"/2 f f ^n^^^^ ^gg^ ^ g^) 

{Jm 



(36) 



+ / d'^-^xVhtl {2QLaa + 12QS + ?>xQ;n) \ • 
JdM ) 

as{Q, L) = -l-(4vr)-("-i)/2 / d-"'x^ii{Q{-2AE + 2AxEx 

+48X^ + 48EX - l2X:aX:a + 12X:aa - QX:aX:aX + 16X^ 
+^XRanan + 1925^ + 96L,,5 + (3 + 10x)i^aa^66 

+ (6 - Ax)La,Lab) + Q;„(965 + 192^2) ^ 24xQ;„„}. (37) 

Since the boundary terms in aQ{l,L) are not known, we have to adopt a 
different strategy to calculate a4(Q,L). The volume part of a^^iQ^L) is already 
known j^lj, so that we have to define the boundary contributions only. First we 
have to write down all possible local boundary invariants of dimension 3 with 
arbitrary coefficients. Boundary invariants are traces over "internal" indices of 
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local polynomials constructed from R, E, Xi ^i L, S and Q and from their 
derivatives. All a,b,c, . . . indices must be contracted in pairs. Note, that the 
normal index n must not be contracted. This reflects specific symmetry of the 
spectral problem in the presence of a boundary which selects direction of the 
normal. Next, we have to use various properties of the heat kernel expansion 
to define these constants. In the particular case Q = If we should restore the 
known result (jMj) . Of course, for a matrix valued Q there are considerably more 
different invariants than in the scalar case since one has to take into account 
non-commutativity of Q with E, x, ^ etc. Therefore, for example, the invariants 
tr (QE) and tr (QxEx) are different although they coincide in the limit Q = If. 
To restrict the number of invariants (and the computational complexity) from 
now on we consider the case 

S = 0, Lab = (38) 

only. 

Since the Riemann tensor does not have internal (spinorial or gauge) indices, 
it commutes with Q and x- Therefore, the particular case Q = If allows to 
restore all curvature dependent boundary terms in a^. 

a^iQ, L) [boundary] = -^(47r)-"/2 f d'^-'xVhtr {Q{l2R,n + 30xi?;n) 

360 Jqm 

+30Q;„X^ + O(i?°)) . (39) 

To control the invariants containing E and VL we use the property [201 ^] 
that the constants appearing in front of all invariants depend on the dimension 
of the manifold only through an overall factor (47?)""^/^. This property makes 
it possible to use low dimensional particular case calculations to define the heat 
kernel coefficients in arbitrary dimension. 



4.2 Particular case calculations 

To define the terms in which depend on x, £" and O it is enough to consider 
the case of the simplest geometry M = M+ x R"~^ with fiat metric. It is im- 
portant that the mass dimension of the boundary integrand in 04 (Q, L) is three. 
Therefore, terms containing both E and Q cannot appear. For this reason, E and 
Q terms can be considered separately. As well, Qat cannot enter the invariants 
since there is no rank two antisymmetric tensor of dimension one. Consequently, 
we may restrict ourselves to the case 

Un = 0. (40) 

To simplify the calculations we also impose 

X = const. (41) 



9 



Then x-.a will be represented by a commutator [cc'ajx], and flan = —dn^a- These 
two invariants are independent on the boundary. 

We shall need a bi-local heat kernel K{x, z\ t) which is defined as a solution 
of the heat equation 

{dt + L) K{x, z;t)=0 (42) 

with the initial condition 

K{x,z;0) = S{x,z) . (43) 
Because of the restrictions and (pOj) the boundary conditions simplify to 

U^K{x,z;t)\9M = 0, dnIl+K{x,z-t)\9M = 0., (44) 
This bi-local kernel is related to the "localised" one (cf. Q) by the equation: 

K{Q,L,t)= [ (Txy/gii {Q{x)K{x,x]t)) . (45) 

J M 

We stress, that K{x, x; t) is a distribution. 

The fundamental solution of "free" heat conduction equation 

{dt-dl)Ko{x,z;t)=0 (46) 

on M = M" is well known: 

Koix'', x", y", y"; t) = (47rt)-"/2exp (- ^°=i^ it ) ' ^^^^ 

From this kernel one can construct a solution of ()46p on M = M+ x M"~^ which 
satisfies the conditions (jl^ at x" = 0: 

K^{x, z; t) = Kq{x\ x", z", t) + xKoix'', x", z\ -z"; t). (48) 

The operators which we consider in this section can be represented in the 
following form: 

L = -5^ - P, (49) 

where P is a first or zeroth order differential operator. 

It is easy to show that the kernel K{x, z\ t) defined by the equation 

K{x,z;t) =K^{x,z;t)+ [ drf dy K^{x,y;t - T)P{y)K{y, z;t) (50) 

Jo J M 

satisfies both the full heat equation (j42|) and the boundary conditions (j44|1 . 
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The equation ()50|) admits a solution^ in terms of power series in P: 



K{x, z; t) = K^{x, z; t) 

+ ^'^P dTp^i... / dri dyp... dyi 

p^;^ Jo Jo Jo Jm J m 

xK^{x, yp] t - Tp)P{yp)K^{yp, ?/p_i; Tp - rp_i) . . . P{yi)K^{yi, z; ri)(51) 

A remarkable feature of (jKTjl is that only a finite number of terms contribute to 
each ak{Q, L) for fixed k. This equation will be used to calculate the heat kernel 
expansion for two particular choices of L (see below). 
We start with a one- dimensional example 

L^ = -dl-Eix), (52) 

so that Pi{y) = E{y). E is an arbitrary matrix valued function on M = IR+. 
Details of the calculation can be found in Appendix El To the first order in E 
we have: 



J.i l>00 / J. \ 

KiQ,L,,t) = -—r / rfxtr QE + -gE/ 
(47r)2 Jo V 6 J 

Q[^iE,n + XE;nX) + \ixE;n + ^;nX)) 

+ 0{t^) . (53) 



(47r)2 

-Q;„(-^(E - xEx) + \ixE + Ex) 



x=0 



The first and the second lines of ()53|) can serve as a consistency check of our 
calculations (cf. (jHHjl . (jHTj) ). The rest of (jK!^ defines uniquely all terms in the 
boundary part of a4{Q,L) containing E subject to the restrictions (jHEI)- Indeed, 
no other invariants of dimension 3 can appear. For example, tr {QE-a) is not 
allowed since it contains a non-contracted tangential index. X;n cannot appear 
since x is defined on the boundary only and, hence, may be differentiated only 
tangent ially. 

In order to define the term in a^^l^Q, L) containing VL and/or tangential deriva- 
tives X we consider a two dimensional example M = Rx]R_|_ and L2 = — (c^^+cij^)^. 
In addition to (pOj) we also suppose daUJa = 0. Clearly, this condition does not 
exclude any relevant invariant. Then P2{y) = 2uja{y^)da + uj1{y^). This time 

^On manifolds without boundary a similar expansion served as a starting point for the 
covariant perturbation theory of ref. j25| . On manifolds with boundary similar perturbation 
series were used [IH] to evaluate dependence of the heat kernel on the boundary conditions. 
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we need the first and second order terms in the expansion ()51|) . All necessary 
technical tools can be found again in Appendix El The result reads: 



M 

t 



dM 

11 3 1 

3U oU 2U Iz 



+ Q 



0(t5). (54) 



Next we collect individual contributions contained in (j^Sl) and (jMj) to 

obtain: 

a4(Q, L) = ^(47r)-"/2| J ^n^^tr {g(60E./ + 60/?^ + 180^^ 
+30^^,^'^^ + 12R,/ + 5R^ - 2R^,R''' + 2R^,p,R^^P'')] 
+ / rf"-^xv^tr{g{30E.„ + 30x^;nX + 90x^;n + 90E.„X 

JdM 

an A :a 

anXi X'-a 

+Q;„(-30E + 30x^X + 90x^ + 90Ex - 

-18X:aX:a + 30X:aa - QX:aX:aX + SO^i?) + SOxQ;/"} }• (55) 



Let us remind that we have imposed a restriction ()38|) on boundary conditions 
and on the extrinsic curvature of the boundary. 



5 Calculation of the anomaly 
5.1 Dimension two 

The remaining part of the calculation is rather simple. One has to substitute 
(jSni) with (j2SD - (IZ3), (El and (HH) in (jHD. The result reads: 

A=~ (^{^-e^^^{V,,+ [A,,A,])-D,M^')^ . (56) 
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Surprisingly, there is no "genuine" boundary contribution here, and this is our 
main result in two dimensions. For constant if it is consistent with an earlier 
calculation The volume term is very well known (see, e.g. [2j). Note, that 
in two dimensions axial vectors can be transformed to vectors, and, therefore, 
the whole anomaly may be generated by just the first term under the integral in 



5.2 Dimension four 

Chiral anomaly in four dimensions is calculated the same way except for that we 
have to use eq. ()55|). Now the anomaly contains two contributions. 



A = A^ + Ab. (57) 



The volume part 



A = ^^//^Xv/^tr^(-120[D,r-,A.] 

+60 [D^A,,V^'']-mD^D^'D,A^ + 12Q{{D^A,,A^},A^} 
+60 {D^A^ A^A''} + 120 A^D.A'^A^ + 30 [[A^, A^A^""] 
+e^upa {-45 1 V^^VP" + 15 i A''" A^"" - 30 1 {V^^AfA" + A^'A^VP") 
-120 1 ApyPA" + 60 i A^A^APA"} - 60 {D„A^)E!"' + 30 {D^A^')R 

-—e,up.R^\eR'''''') (58) 

is known^. 

The boundary part 

•^^ = T^?77^ / d^xVhtrluie''''' {A,,ip}DaA, 
180 (27r)^ Jqm ^ 

+24{(^, An{Aa, A^} - 60 [A^ if] (Ka - [An, A]) 

+60{Dnip)D^A^) (59) 



is new. It has been derived under two restrictions (j38|) . Note, that in the present 
context, the first condition (5* = 0) actually follows from the second one [Lab = 0) 
due to ((221) • We shall analyse physical consequences of in a future publica- 
tion. 



^The flat space part of (|58|) can be found e.g. in The term with R-k R follows from the 
local index theorem and has a very long history |27j . We are not aware of any works which 
considered chiral anomaly for non-abelian axial vector field in curved space, so that the terms 
{DA)R have a chance to be new. There are certain similarities between H58|l and chiral anomaly 
in the Riemann-Cartan space |28| . 
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6 Conclusions 



We have calculated boundary contributions to cliiral anomaly for local (bag) 
boundary conditions in two and four dimensions (in four dimensions we have 
supposed that the boundary is totally geodesic, Lab = 0). As a by-product 
we obtained explicit expressions for several heat kernel coefficients with mixed 
boundary conditions and with a matrix-valued smearing function. These heat 
kernel expressions are rather universal. By choosing a bit different expressions 
for Q and for the fields - (j?fj) one can easily extend our results to other 
anomaly-like expression relevant for hadron physics (see, e.g., [2^]) and, probably, 
even to supersymmetry [HU] . 

We have found no specific boundary contributions to the anomaly in two 
dimensions. In four dimensions there are boundary terms in the anomaly, which 
must have important physical consequences both in hadron physics and in the 
standard model. 

Our present paper is the first one which treats local chiral anomaly in the pres- 
ence of boundaries in generic background vector and axial vector fields. Therefore, 
our results may be improved or extended in many directions. 

It is clear that in order to bring our results closer to physical applications 
we have to lift the restriction Lab = in four dimensions (which excludes, for 
example, spherical boundaries in fiat space). This is just a technical problem 
which can be solved by the same methods as presented above. Another problem 
is to extend our results beyond the linear order in the chiral transformation 
parameter. Such an extension requires chirally transformed boundary conditions 
fll7|l . We are going to address these two problems in the near future. 

Brane-world and domain wall configurations lead to interactions confined at a 
singular surface. Mathematically such interactions are described at the one-loop 
level by some "gluing conditions" which relate boundary values of the functions 
and their normal derivatives on two sides of the singular surface. Such cases 
may be treated by the same methods as presented in this paper. Moreover, a 
lot of important information on the heat kernel expansion for gluing conditions 
is contained in the heat kernel expansion for the boundary conditions case (see, 
e.g., ED). 
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A Heat kernel expansion with mixed boundary 
conditions 



Here we give explicit expressions for the heat kernel coefficients a^, A; = 0, 1, 2, 3, 4 
for an operator of Laplace type subject to mixed boundary conditions [20] (see 
also j21j for minor corrections in 04), / is a scalar function. 



ao{f,L) = (47r)-"/^ / (Tx^iiif) (60) 

J M 

a,{f, L) = i(47r)-("-i)/2 f d^-^x^^v (x/). (61) 

^ JdM 

a^if, L) = ^(47r)-"/2 rx^tr {6fE + fR) 

+ [ r-'xVhtTi2fLaa + Sxf-n + l2fS)] . (62) 
JdM J 

asif, L) = -i^(47r)-("-i)/2 / d"-^xv^tr {/(96xi? + l^xR 

+ SfxRanan + (130+ - 7U_)LaaLtb + (2n+ + 10n_)L,f,L,, 
+ 96SLaa + 192^2 - 12X:aX:a) + /;n((6n+ + 30n_)L,, 

+96^) + 24x/;nn}. (63) 
a,if,L) = :^{4n)-^/\ [ d^x^ti {f{6QE,^'^ + GORE + 

+30^^,^^^" + 12R./ + bR" - 2i?^,i?^" + 2i?^,p,i?'^"'"^)} 
+ / rf"~W/Itr{/{(240n+ - 120n_)E;„ 

JdM 

+ (42n+ - 18n_)i?.„ + 24LaaM + l20ELaa 

+ 20RLaa + 4:RananLbb — ^'^RanhnLab + ^RahcbLac 

+ ^{(280n+ + AOYiJ)LaaLbbLcc + (168n+ 
-2UIl.)LabLabL^ + (224n+ + 320n_)L„bLfe,L,J 
+720^E + 120SR + 144^L„„Lfe;, + ASSLabLab 

+A%OS''Laa + 4805=^ + 1205:,, + mXXMan ' l^X.aX-.aUb 
-2AX:aX:bLab ~ l20X:aX:aS} + /;„(180x^ + 30^/? 

+ ^{(84n+ - 180U_)LaaLbb + (84n+ + QOU_)LabLab} 

+ 72SLaa + 240^2 - 18X:aX:a) + f;nn{2ALaa + 1205) 

+30x/;/"}}. (64) 

The coefficient can be found in [3^. Here we present only those terms 
which enter ()32|) with A; = 3. Namely, we put / = 1 and neglect all terms which 



15 



do not contain E. 

a,{l, L) = -i-(47r)-(™-i)/2 / d^-'xVhti (360x^;„„ + 1440E;„5 
5760 Jqm 

+72QxE^ + l80x:aaE + 2A0xRE - UOxRnnE + 2880ES^ 

+ (270 - l80x)L,aE.n + lUOLaaSE + (45 + l50x)LaaLbbE 

+ (90 - QOx)LabLabE - 180(^2 - xExE) - ISOx-.aX-.aE 

-90xX:aX:aE)+OiE'). (65) 



B Particular case calculation: details 

Let us consider in detail calculations of the expansion ()53|) for the operator ()52|) . 
First we use the fl51|) to represent the first order terms in E of K{Q, Li] t) as 

t OO CO 

1- j dr j dx j dy-^^^===tT !^Q{x) (^e"^ + xe"^ 



X E{y) e + 4{t^ ^ (66) 



Next we integrate over r with the help of the relation 



',/M~f:><'M^&^^.^j\A±Mi (67) 



The parameters a and 6 are equal either to {x — y)/2 or to (x + y)/2. Note, that 
"reflected" terms depending on (x + y) are always multiplied by x- 

It is convenient to consider 3 different types of contributions separately. 

1. Terms without x- 

In this case a = b = {x — y)/2. Let us change the variables 

X = k + rVi, y = k for X > y 

y = k + rVi, x = k for x < y. (68) 

In both cases k,r & [0, +oo[. Then we use a small t expansion: 

r..erfcW/(.v/F) = /p + n2M + + (69) 

Jo VTT 4 b^y7^ 

which is valid for any smooth function which decays sufficiently fast at 
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infinity. We obtain: 

-OO /• + 00 



+00 .+00 (g^^ ^ Vtr)Eik) + g(A;)^(A; + v^r)) 

/ ak dr eric r tr 

Jo Jo 4 

dytrl Q{y)E{y) + -Q{y)E"{y) - -trQ(0)E(0) 



7o V 6^'^' '^7 16 

+ ^tr (^g(O)i? '(0) - '(O)E(O)) + Oie) . (70) 

2. Terms with two x- this case a = b = {x + y)/2. We change the variables: 

x = rv^cos0, y = r\/i sin (f) , (71) 

so that r G [0, oo[, e [0, 7r/2]. We use (jHIj) and then integrate over r and 
(j) assuming t is small . The resulting asymptotic expansion reads 

3 

+^tT[^Q'{0)xE{0)x+^QmE'm)+O{t') (72) 



3. Terms with single x- We have a = {x — y)/2, b = (x + y)/2. We use the 
variables: 

X = r-\/t(cos + sin 0), y = r\/i cos (j) for x > y; 

y = rVi{cos(f) + sm(j)), x = r^/i cos (f) ioix<y. (73) 

Acting as above we obtain: 

=° , , . / |x + y| + |x-y| V g(x)(xi^(i/) + E(y)x) 

dx / ay eric ^ tr ^ = 

Jo \ 2Vi J 4 

UtQ{0){xE{0)+E{0)x) + 

3 

^tr (ig(0) (xi? '(0) + E '(0)x) + ig '(0) {xE{0) + E{0)x) ) + 0{t') 

(74) 

The sum of dZOI), (IZ21) and (HH) yields 
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Calculations of the uj terms can be carried out in a similar way. The following 
integral is useful: 



* rv exp(-^)exp(--^)exp(-^) 

ctp / dr — 
Jo 



\Jt - V\/V - 

-jb\_ 

Vt 



-27r'^/^-(|a| + |6| + |c|)-erfc ^— ^ — — +27rv^-expf ^ ' 



(75) 
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